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We describe the design for a scalable, solid-state quantum-information-processing architecture 
based on the integration of GHz-frequency nanomechanical resonators with Josephson tunnel junc- 
tions, which has the potential for demonstrating a variety of single- and multi-qubit operations 
critical to quantum computation. The computational qubits are eigenstates of large-area, current- 
biased Josephson junctions, manipulated and measured using strobed external circuitry. Two or 
more of these phase qubits are capacitively coupled to a high-quality-factor piezoelectric nanoelec- 
tromechanical disk resonator, which forms the backbone of our architecture, and which enables 
coherent coupling of the qubits. The integrated system is analogous to one or more few-level atoms 
(the Josephson junction qubits) in an electromagnetic cavity (the nanomechanical resonator). How- 
ever, unlike existing approaches using atoms in electromagnetic cavities, here we can individually 
tune the level spacing of the "atoms" and control their "electromagnetic" interaction strength. We 
show theoretically that quantum states prepared in a Josephson junction can be passed to the 
nanomechanical resonator and stored there, and then can be passed back to the original junction or 
transferred to another with high fidelity. The resonator can also be used to produce maximally en- 
tangled Bell states between a pair of Josephson junctions. Many such junction-resonator complexes 
can assembled in a hub-and-spoke layout, resulting in a large-scale quantum circuit. Our proposed 
architecture combines desirable features of both solid-state and cavity quantum electrodynamics 
approaches, and could make quantum information processing possible in a scalable, solid-state en- 
vironment. 

PACS numbers: 03.67.Lx, 85.25.Cp, 85.85.+j 



I. INTRODUCTION 

The lack of a large collection of physical qubit ele- 
ments, having both sufficiently long quantum-coherence 
lifetimes and the means for producing and controlling 
their entanglement, remains the principal roadblock to 
building a large-scale quantum computer. Supercon- 
ducting devices have been understood for several years 
to be natural candidates for quantum computation, 
given that they exhibit robust macroscopic quantum 
behavior.^- Demonstrations of long-lived Rabi oscilla- 
tions in current-biased Josephson tunnel junctions, ^'^ 
and of both Rabi oscillations and Ramsey fringes in 
a Cooper-pair box,**'^'^ have generated significant new 
interest in the potential for superconductor-based 
quantum computation^*^ Several additional experimen- 
tal accomplishments have followed|^i^°i^^i^^i^^i^'^i^^i^^ 
including the impressive demonstration of controUed- 
NOT logic with charge qubits, and a large body of 
theoretical work is beginning to address these and related 
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Coherence times r^p up to 5 ^s have been reported in the 
current-biased deviceSfSi with corresponding quantum- 
coherent quality factors Q^p = AE/h of the order 
of 10^, indicating that these systems should be able to 
perform many logical operations during the available 
coherence lifetime Here AE is the qubit energy-level 
separation, which was 68 fieV in the experiment of 
Ref. [2. 

In this paper, we expand on our earlier proposal sug- 



gesting that GHz-frequency nanoelectromechanical res- 
onators can be used to coherently couple two or more 
current-biased Josephson junction (JJ) devices together 
to make a flexible and scalable solid-state quantum- 
information-processing architecture^ The computa- 
tional qubits are taken to be the energy eigenstates of 
the JJs, which are to be individually prepared, controlled, 
and measured using the external circuitry developed by 
Martinis et These superconducting phase qubits are 
capacitively coupled to a high- quality- fact or piezoelectric 
dilatational disk resonator, cooled on a dilution refriger- 
ator to the quantum limit, which forms the backbone of 
our architecture. We shall show that the integrated sys- 
tem is analogous to one or more few-level atoms (the JJs) 
in an electromagnetic cavity (the resonator). However, 
here we can individually tune the energy level spacing of 
each "atom" , and control the "electromagnetic" interac- 
tion strength. This analogy makes it clear that our design 
is sufficiently flexible to be able to carry out essentially 
any operation that can be done using other architectures, 
,37.3BmyMad4dijEt4feibtH:e6iax:«>ug to.6(D herence. Many of our re- 
sults will apply to other architectures that are similar to 
atoms in a cavity. 

Several investigators have proposed the use of 
LC resonators superconducting 
cavities or other types of oscillators^SiSii^^ to 
couple JJs together. We note that although harmonic os- 
cillators are ineffective as computational qubits, because 
the lowest pair of levels cannot be frequency selected by 
an external driving force, they are quite desirable as bus 
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qubits or coupling elements. Early on, Shnirman et alm^ 
suggested an architecture consisting of several supercon- 
ducting charge qubits in parallel with an inductor. The 
J Js are themselves out of resonance with each other and 
hence weakly coupled, and the resulting LC resonator 
(the capacitance coming from the junction geometry) is 
also used well below its resonant frequency. An interest- 
ing modification of this design couples the small island to 
the external circuit through a pair of parallel JJs, which 
allows the Josephson coupling energy to be varied, us- 
ing an external magnetic field. To date, however, the 
only coupled superconducting qubits demonstrated ex- 
perimentally have been the capacitively coupled charge 
qubits of Peshkin et al^ and Yamamoto et al.^ and the 
capacitively coupled phase qubits of Berkley et alm^ 

Resonator-based coupling schemes, such as the one 
proposed here, have the advantage of additional func- 
tionality resulting from the ability to tune the qubits rel- 
ative to the resonator frequency, as well as to each other. 
We shall show that by tuning the JJs in and out of res- 
onance with the nanomechanical resonator, qubit states 
prepared in a junction can be passed to the resonator and 
stored there, and can later be passed back to the original 
junction or transferred to another JJ with high fidelity. 
The resonator can also be used to produce controlled 
entangled states between a pair of JJs. Alternatively, 
when both qubits are detuned from the resonator, the 
resonator produces a weak (higher order) "dispersive" 
qubit coupling similar to that of a capacitor. The use of 
mechanical resonators to mediate multi-qubit operations 
in JJ-based quantum information processors has not (to 
the best of our knowledge) been considered previously, 
but our proposal builds on the interesting recent theo- 
retical work by Armour et alm^^ and Irish et alm^ on 
the entanglement of a nanoelectromechanical resonator 
with a single Cooper-pair box. In fact, there is currently 
a big effort to push a variety of nanomechanical systems 
to the quantum limit iSiS^iS 

In the next section we recall the basic properties of 
large-area, current-biased JJs. In Sec. Illll we discuss 
our proposed architecture, and construct a simple model 
Hamiltonian to describe it. State-preparation and read- 
out have been described elsewhere and will only be dis- 
cussed briefly. The properties of the nanomechanical res- 
onator are also described here in detail. In the remainder 
of the paper we discuss a variety of elementary single- and 
multi-qubit operations central to quantum computation: 
In Secl^we show how a qubit state prepared in a J J can 
be passed to the nanomechanical resonator, stored there 
coherently, and later passed back to the original junc- 
tion or transferred to another JJ. Two-junction entan- 
glement, mediated by the resonator, is studied in Sec.lVI 
In Sec. lVII we show how our architecture can be extended 
to make a large-scale quantum circuit. Our conclusions 
are given in Sec. IVIII Several immediate extensions of 
the present work, including the development of protocols 
for universal two-qubit quantum logic, are currently in 
progress and will be discussed in future publications. 



II. THE CURRENT-BIASED JOSEPHSON 
JUNCTION 

Our architecture relies on the use of large-area JJs, bi- 
ased with a current lb, which can be quasi-static or have 
oscillatory components. The junctions have a large ca- 
pacitance C (typically 1 to 50 pF) and critical current 
Iq (in the 10 to 150 /iA range) so that the largest rele- 
vant energy scale in the system is the Josephson coupling 
energy 



E, = 



2e 



(1) 



where e is the magnitude of the electron charge. In con- 
trast, the Cooper-pair charging energy 



2C 



(2) 



is small compared with Ej, and is also usually smaller 
than the thermal energy ksT. For example, 

Ej = 2.05 meV x /q [/iA] and E, - ^^J^ , (3) 

C[pt\ 

where /o[AiA] and C[pF] are the critical current and junc- 
tion capacitance in microamperes and picofarads, respec- 
tively. 



A. Semiclassical junction dynamics 

The dynamics of a real JJ can be understood as fol- 
lowing from the equivalent circuit model shown in the 
inset to Fig. ^ known as the resistively and capacitively 
shunted junction modelj^SiS where the "ideal" Joseph- 
son element controls the superconducting component Is 
of the total electrical current Is + In in accordance with 
the well-known Josephson equations 



Iq sin 6 and 



dt 



2eV 



(4) 



Here Is is the supcrcurrent flowing through the ideal 
Josephson junction element, S is the difference between 
the phases of the (spatially uniform) superconducting or- 
der parameters on each side of the junction, and V is the 
voltage across the junction. /„ is the normal compo- 
nent of the current. The resistance R accounts for finite- 
temperature quasiparticle tunneling as well as electron 
tunneling in the flnite-voltage state. Equating the sum 
of the currents flowing through the capacitor, ideal junc- 
tion, and resistor, to /b, leads to 



2Ec dt^ ^ Ae'^R dt 



dS 



sin (5 



= 0, (5) 



where s — Ib/Io is the dimensionless bias current. By 
rewriting Eq. ^ in the equivalent form 



M 



df^ 



dU 



V 



dS 

dt' 



(6) 
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it can be interpreted as the equation of motion for a 
particle of "mass" 



B. Quantizing the low-energy junction dynamics: 
The phase qubit 



M = 



2E, 



(7) 



(8) 



moving in an effective potential 

U{S) = -Ej{cosS + sS), 

and in the presence of velocity-dependent dissipation 
characterized by a friction coefficient 77 = H^/Ae^R. Note 
that M actually has dimensions of mass x length^. 

The potential U{6), which resembles a tilted wash- 
board, is shown in the main panel of Fig. ^ for a dimen- 
sionless bias current of s = 0.1. The zero- voltage state of 
the junction corresponds to the particle or phase variable 
being trapped in one of the metastable minima present 
when s < 1, and the finite- voltage state corresponds to 
the phase variable running down the washboard poten- 
tial. In what follows we will assume < s < 1, and with- 
out loss of generality we can also assume that < S < 2tt. 

The potential U{5) reaches its minimum and maximum 
values in the domain < 5 < 27r at (Jmin — arcsins and 

<^max = TT - arCSiu S. The depth AC/ = C/((5max) - U ((5,nin) 

of the potential well is 



AC/ = 2Ej vl — — s arccos s , 



which vanishes as 



AC/ 



4%/2 



Ei{l-s) 



3/2 



(9) 



(10) 



in the s ^ 1" limit. 

The curvature U"{6) at the minimum of the potential 
is used to define the junction's plasma frequency, 



[/"(-Jn 



M 



^po(l 



(11) 



which is the frequency of small oscillations of 6 about 
(5min- Here M is the effective mass defined in Eq. {TJ, 
and 



2e/o V2KE1 



hC 



(12) 



is the plasma frequency at zero bias. 

The dependence of the barrier height and plasma fre- 
quency on bias current are plotted in Fig. |21 For junc- 
tions appropriate for quantum computation, ti;po/27r is 
typically in the range of 5 to 50 GHz. The barrier height 
during state preparation and readout is usually adjusted 
so that AU/hu!p is between 3 and 5, but, as we shall dis- 
cuss below, is it advantageous to keep s smaller during 
actual quantum computation. 

The effect of dissipation, caused in the resistively and 
capacitively shunted junction model by the resistance R, 
can be characterized by the number of oscillations at 
the plasma frequency during the relaxation time RC, or 
LUpRC. In what follows we will assume that dynamics is 
highly underdamped, with ujpRC 3> 1. 



When the thermal energy fcsT and energy decay width 
h/RC are both smaller than tuvp, quantum fluctua- 
tions of 6 become important, and the J J has to be 
treated quantum mechanically. This limit was stud- 
ied in the 1980's as an example of a single macro- 
scopic degree of freedom — the difference between phases 
of order parameters — that nonetheless behaves quantum 
mechanicallyeSSiSMMi This is also the regime of current 
interest for applications to quantum computing. 

When dissipation is absent, the low-energy dynam- 
ics can be quantized by introducing a Lagrangian Lj = 
i/V/j^ — U and canonical momentum P = dL/d6 — MS 
associated with the 77 = limit of Eq. © . According to 
the Josephson equations, P is proportional to the charge 
Q or to the number of Cooper pairs Q/2e on the capaci- 
tor according to P = hQ/2e. The classical Hamiltonian is 
P^/2M + U. To quantize the system, we let P 
so that [S, P] — ih. Then the quantized Hamiltonian is 



(13) 



and the dynamics is governed by the Schrodinger equa- 
tion ihdtip — Hjip. Because U depends on s, which itself 
depends on t, Hj is generally time-dependent. 

Naively, the stationary states and energies of the 
JJ with fixed s follow straightforwardly from the one- 
dimensional eigenvalue problem 



Hj ^miS) ^ em1pm{S), TO = 0,1,2, 



(14) 



However, a careful analysisSS shows that in the presence 
of any finite ohmic dissipation (finite rj), quantum co- 
herence between the different wells in U{S) is destroyed. 
This, in fact, justifies the use of the washboard potential 
in the first place: Strictly speaking, (5 is a periodic vari- 
able, with S physically equivalent to i5 -|- 2tt. In what fol- 
lows, we will work with stationary states associated with 
a single potential minimum (in the domain < i5 < 27r). 
It is these stationary states that are of interest to quan- 
tum computation. 

When s = 0, the junction contains many (of order 
^Ej/Ec) bound states, the lowest of which are like that 
of a harmonic oscillator with level spacing hiOpQ . The uni- 
form spacing of the low-lying levels makes them difficult 
to address individually with a classical external driving 
force. Therefore, state preparation is carried out with s 
just below unity, in which case there are only a few qua- 
sibound states |0), |1), |2), . . . present, and the effective 
potential U{6) becomes anharmonic and approximately 
cubic, as illustrated in Fig. |31 The remarkable 1985 
spectroscopic observation^ of these quantized states pro- 
vided the first clear evidence for the quantum behavior 
of the macroscopic phase-difference variable 6. 

The lowest two eigenstates, |0) and |1), define a phase 
qubit. As stated, in the s ^ 1 limit the potential is 
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anharmonic, and the qubit level spacing 

AE = ei- eo (15) 

is somewhat smaller than fkOp, where cOp is the s- 
dependent plasma frequency. 

The qubit state is also usually measured with s just be- 
low unity: In the absence of thermal or quantum fluctu- 
ations, switching to the finite-voltage state occurs when 
the bias current exceeds Iq. However, in a real junc- 
tion, the finite-voltage state will occur before Ih reaches 
/q, either because of thermal activation over the barrier 
or by quantum tunneling through it. Once the phase 
variable escapes into the continuum, it runs down the 
corrugated potential, and a voltage V of approximately 
2Asc/e develops across the junction, where Age is the 
superconducting energy gap (Age ~ 180 ^eV for Al junc- 
tions) . The supercurrent component then oscillates with 
angular frequency 2eV/h — the AC Josephson effect. The 
thermal activation regime has been explored in detail, 
for various limits of dissipationSSM^ For fixed current 
bias, the thermal activation rate falls exponentially with 
inverse temperature, until the dominant escape mech- 
anism becomes quantum tunnelingi^2iS24SSiSLS& At tem- 
peratures low enough so that quantum tunneling domi- 
nates thermal activation, the qubit state can be observed 
by measuring the tunneling rate, which is strongly state- 
dependent. State preparation and readout are discussed 
further in Sec. IIII Al 

The barrier height AU and the energy splitting AE 
(through its dependence on ujp) are both strong functions 
of the bias current s. The ability to tune the plasma 
frequency is one of the current-biased Josephson junc- 
tion's great strengths and weaknesses. It enables the 
qubit level-spacing AE to be tuned adiabatically into 
resonance with another qubit or, as in our approach, 
with a resonator, but it also makes the circuit sensitive 
to bias-current noise, as characterized by the non-zero 
derivative dup/ds. Fluctuations in s will generate noise 
and hence decoherence in the J J. Although current meth- 
ods of state preparation and measurement require s very 
close to unity (typically near 0.99), where dojp/ds is un- 
fortunately large, the information-processing operations 
we describe below do not. In our simulations, we find it 
convenient to work with s below 0.90. 

The energies e™ of the lowest four JJ states of the de- 
vice used in Ref. for a range of bias currents, are 
given in Table in units of hiUp. We calculate these ener- 
gies numerically by diagonalizing the Hamiltonian Hj of 
Eq. (|13|l in a basis of harmonic oscillator eigenfunctions 

= (2™m!V^4)-'/' e-«'/2 Hm{0 (16) 

that are constructed by making a quadratic approxima- 
tion 

U{6) « U{5ra^r.) + \ U"{6„,in){S - ^min)' (17) 

to U{S) about its minimum at dmin = arcsins. The 
Hm, m = 0, 1, 2, . . . , are Hermite polynomials, and ^ = 



TABLE I: Energies of low-lying eigenstates as a func- 
tion of dimensionless bias current s, for the JJ investigated 
in Ref. 3j, with parameters Jo = 21 A {Ej — 43.05 meV) 
and C — 6pF {Ec — 53.33 neV). Energies below are given 
in units of hcup and are measured relative to i7(Jmin). All 
dissipation and decoherence effects are neglected. The first 
column, labeled by m + i, gives the energies of the corre- 
sponding harmonic oscillator eigenfunctions, which are found 
to be extremely accurate for small s. 



junction state |m) 




s = 0.50 


s = 0.70 


s = 0.90 


m = 


0.500 


0.500 


0.500 


0.500 


m — 1 


1.500 


1.500 


1.499 


1.497 


m = 2 


2.500 


2.499 


2.498 


2.492 


m — 3 


3.500 


3.498 


3.496 


3.485 



TABLE II: Dipole moments Xmm' between pairs of low-lying 
JJ eigenstates states for bias s — 0.90. The entries with dots 
follow from symmetry. Junction parameters are the same as 
in Table B 



{m\S\m') 


m' = m' = 1 


m' = 2 


m' = 3 


m = 


1.12 3.46x10"^ 


-5.86x10"* 


7.09x10"'^ 


m = 1 


1.12 


4.89x10"^ 


-1.02x10"^ 


m = 2 




1.13 


6.00x10"^ 


m = 3 






1.13 



{S — dmin)/£s is a Centered and scaled phase variable, with 




giving the characteristic width in 6 of these eigenfunc- 
tions. We find rapid convergence to the values reported 
in Table ^ as the number of harmonic oscillator basis 
states is increased to include all basis states with ener- 
gies less than U (5max) • 

Dipole-moment matrix elements 

Xmm' = {m\S\m'), (19) 

which will also be used below, are calculated at bias 
s = 0.90 for the junction used in Ref. Q, using this same 
method. The results are given in Table Hll All basis func- 
tions with energies less than J7(5max) are included, and 
the oscillator-strength sum rules (adapted for this Hamil- 
tionian) are satisfied to better than 99.999%. Because 
the eigenfunctions are real, the matrix Xmm' is symmet- 
ric, and with an appropriate choice of overall signs of the 
eigenfunctions, the first band of off-diagonal matrix ele- 
ments can be made positive. The diagonal elements are 
also positive here, a consequence of our restriction to the 
domain < (5 < 27r. 

The diagonal elements Xmm are very close to Smin, re- 
gardless of TO. In the s = 0.90 case considered in Table Ull 
i5inin is about 1.120. The values of off-diagonal elements 
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of the form Xm,m±i can be understood by noting that for 
harmonic oscillator states, which in this case are close to 
the exact eigenf unctions, 

JdS <P^{S)S^^+,iS) = ^1^^ (20) 

with £s = 4.883 X 10^^. The remaining off-diagonal el- 
ements, which result from the small mixing of the har- 
monic oscillator states, are smaller than these by at least 
an order of magnitude. 



architecture consisting of a network of current-biased 
Josephson junctions coupled to nanoelectromechanical 
resonators. We will first consider a single nanomechanical 
resonator coupled to one or two JJ qubits; the extension 
to larger systems will be considered below in Sec. IVII as 
well as in future work. 



III. ARCHITECTURE AND MODEL 
HAMILTONIAN 

We turn now to the main focus of our paper, the de- 
scription of a solid-state quantum-information-processing 



The complete circuit diagram for the two-JJ circuit 
is shown in Fig. 01 The two central crossed boxes are 
the JJs to be used as phase qubits, and they include the 
parallel capacitance and resistance shown in the inset to 
Fig. n] The disk-shaped element in the center of the 
figure is the nanomechanical resonator, consisting of a 
single-crystal piezoelectric disk sandwiched between two 
metal electrodes. Applying a voltage across this element 
produces an electric field between the plates, and through 
the piezoelectric response, a strain in the crystal. Con- 
versely, strain in the resonator produces a charge on the 
electrodes, whose rate of change contributes to the cur- 
rent flowing through a J J. 



A. Single-qubit state preparation, manipulation, 
and readout 

Two of the most critical factors in the design of a suc- 
cessful JJ-based quantum information processor are high- 
impedance bias and high-fidelity readout circuits that do 
not disturb the qubit during computation. This is cur- 
rently a subject of active experimental investigation, and 
for concreteness we will assume the bias circuit design de- 
veloped recently by Martinis et al.f but we will leave the 
the readout circuitry unspecified. Our architecture can 
be adapted to improved readout schemes as they become 
available. 

State preparation and readout are performed with /b 
just below /q, where U{S) is anharmonic and shallow. 
The anharmonicity allows preparation from a harmoni- 
cally varying bias current, which is tuned to couple to 
only the lowest two states. The |0) state is prepared by 
waiting for any excited component to decay. The state 
or a superposition a\0) + is prepared by adding 
radiofrequency (RF) components of magnitudes I^^ and 



I^f to the DC bias current, in the formai 

Ih{t) = Idc + /rf cos{ujrit) + /rf sin(tjrfi), (21) 

with /(jc and I'^^'^ all varying adiabatically (slow compared 
with the frequency AE/H). When oj^i is nearly resonant 
with AE/h, the qubit will undergo Rabi oscillations, al- 
lowing the preparation of arbitrary linear combinations 
of |0) and The associated Rabi frequency 




depends on both Srf = Irt/Io and the detuning. All states 
on the Bloch sphere may be prepared in this manner 41 
Readout of a JJ state a\Q) + (3\1) is performed by then 
tuning into resonance with (e2 — ei)/^-, thereby excit- 
ing the qubit component in the |1) state up to |2), out 
of which it quickly tunnels, thereby resulting in a mea- 
surement of |/3|. Martinis et al^ have established that 
a single-shot readout of the JJ states |0) and |1) can be 
performed with 99% and 85% accuracy, respectively. 



B. Nanomechanical resonator 

The second important element in our design is the 
use of piezoelectric nanoelectromechanical disk res- 
onators, with dilatational-mode frequencies u)q/2tt in the 
1 to 50 GHz range. Piezoelectric dilatational resonators 
with frequencies in this range, and quality factors Q = 
loqt of the order of 10"^ at room temperature, have been 
fabricated from sputtered AlNiS2i2& Here t is the energy 
damping time. The radius of the disk is denoted by R, 
and b is its thickness. In Ref. [sj we presented resonance 
data down to 4.2 K for a 1.8 GHz AIN resonator. The 
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observed low-temperature Q of 3500 corresponds to an 
energy lifetime r of more than 300 ns, already sufficient 
for most of the operations described below. This is to 
be contrasted with the previous state-of-the-art 1 GHz 
SiC cantilever beam resonator demonstrated in 2003/^ 
which has a Q nearly an order of magnitude smaller at 
the same temperature. The unprecedented performance 
of our resonator is a consequence of the use of AIN, which 
is an intrinsically high Q material)^^ and the use of the 
dilatational vibrational mode. 

The dilatational mode of interest is an approximately 
uniform oscillation of the thickness of the disk, which 
produces a nearly uniform electric field in a direction 
perpendicular to the disk. For a disk with large aspect 
ratio R/h, the dilatational mode frequency is 

cjQ = TTv/b, (23) 

where u is a piezoelectrically enhanced sound speed to 
be defined below. Although the dilatational mode is not 
necessarily the fundamental mode of the resonator, we 
can couple to it by frequency selection, carefully avoid- 
ing the other low- frequency modes. The frequency in 
Eq. (|23|l is that of the fundamental vibrational mode of 
a one-dimensional elastic string with free ends. For sim- 
plicity, we will assume that the dilatational mode fre- 
quency given by Eq. (|23|) holds even if the aspect ratio 
R/h IS not large. 

Quantum mechanically, each vibrational mode n of 
such a resonator, having angular frequency w„, is equiv- 
alent to a harmonic oscillator with energy level spacing 
fujjn. For sufficiently high frequency and low temper- 
ature, the mode can be cooled to its quantum ground 
state: For example, if ujqI1-k = 15 GHz, then hio^jk^ 
is about 720 mK. If cooled on a dilution refrigerator to 
100 mK, the probability 

pi = 2 sinh(2|^) e-3'"«^o/2feBT (24) 

of thermally occupying the first excited (one-phonon) 
state, thereby producing a mixed state instead of the 
desired pure phonon ground state, is smaller than 10"^. 
The mean number nTg,(hwQ) of phonons present in the di- 
latational mode at 100 mK, or "excitation level" of the 
corresponding harmonic oscillator, is also less than 10~^. 
Here nB(e) is the Bose distribution function. 

In the simulations below we will assume a nanomechan- 
ical disk resonator with the parameters given in Table lTTll 
The thickness h is chosen to give a dilatational mode fre- 
quency Wo/27r of 15 GHz. This frequency is convenient for 
simulation because, when coupled to a JJ with param- 
eters corresponding to that of Ref. ^ , the bias current 

s* = .^l-(c^oKo)4 (25) 

required to tune the qubit level spacing Ai? into reso- 
nance with hiOQ is small enough so that the JJ eigenfunc- 
tions can be taken to be harmonic oscillator states. The 



TABLE III: Parameters characterizing the piezoelectric res- 
onator simulated in this paper. 



piezoelectric material 


AIN 


mass density p 


3.26 gcm~^ 


dielectric constant £33 / eo 


10.7 


elastic stiffness C33 


395 GPa 


piezoelectric modulus 633 


1.46 Cm"^ 


piezoelectric efficiency ^ == 633/633033 


0.057 


enhanced stiffness £33 = (1 -I- 7)033 


418 GPa 


sound velocity v = C33 /p 


11.3 kms"^ 


disk radius R 


0.230 pm 


disk thickness b 


377 nm 


dilatational frequency ujof^n 


15 GHz 


frequency in Kelvin hcjo/kB 


720 mK 


resonator capacitance dcs 


0.042 £F 



resonator radius R listed in Tabic ITTll is chosen to make 
the junction-resonator interaction strength g, to be de- 
fined below, 1% of hujQ, although we will also briefly con- 
sider larger resonators with larger interaction strengths. 
As stated above, the parameters listed in Table IIIII as- 
sume that Eq. (|23|l is valid. The AIN physical constants 
were obtained from the review by AmbacherJ'^ 

We turn now to a calculation of the dilatational mode 
of the piezoelectric disk, assuming & <C i?. The disk lies in 
the xy plane. In the b <^ R limit, the elastic displacement 
field u(r, t) for the dilatational mode is directed in the z 
direction, and the z component is itself only dependent 
on z and t. Edge effects are assumed to be negligible. The 
vibrational dynamics for this mode and its harmonics is 
therefore effectively one-dimensional. 

Let u denote the z component of the displacement 
field. To construct the equation of motion for u(z,t), 
we write the basic electromechanical equations of piezo- 
electric media^** in the form 

= —D, - d,u, (26) 
£33 

and 

T,,^-h3sD,+d33d,u. (27) 

Here and Dz are the z components of the electric 
E and D fields, and is the stress tensor. €33 is 
the relevant element of the static dielectric tensor, and 
/i33 = 633/633, with 633 the piezoelectric modulus. Fi- 
nally, C33 = (1 -|- 7)033 is a piezoelectrically enhanced 
elastic modulus, with C33 denoting the appropriate ele- 
ment of the elastic tensor, and 



633633 

is a dimensionless quantity called the piezoelectric effi- 
ciency. The values of these material parameters for the 
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case of AIN are summarized in Table HTll Eq. 1)26(1 deter- 
mines the relation between the electric field and strain 
inside the resonator, and Eq. H27|l determines the stress- 
strain relationship, as modified by the electric field. 

Electrically, the boundary conditions are that there is 
a charge per unit area a on the top electrode of a parallel- 
plate capacitor enclosing the resonator, and —a on the 
lower electrode. Then, in the interior of the piezoelectric, 
Dz is uniform, with the value 



returning to the semiclassical description of the JJ re- 
viewed in Sec. Ill AI including the resonator in the equiv- 
alent circuit, and then requantizing the coupled system. 

Our first objective is to derive an equation for /res, 
the resonator's contribution to the electrical current seen 
by the JJ. /res is equal to g, where q is the charge on 
the resonator's top (ungrounded) electrode produced by 
voltage fluctuations across and strain fluctuations inside 
the resonator. Integrating Eq. H26() gives the voltage 



(29) 



Mechanically, the faces of the resonator are assumed to be 
stress free. We note from Eq. H27|l that when cr 7^ 0, this 
stress-free condition requires a fixed strain of ~h33a/c:i3 
on the upper and lower surfaces of the disk. Note that 
these boundary conditions are generally time-dependent, 
because a usually is. 

The resonator has thickness b and occupies the region 
< z < 6. From the mechanical equation of motion 
p&fui = —djTij we obtain 



{dt - v^dl)u = 0, with V = VeiWp- (30) 

The sound velocity in the z direction is slightly enhanced 
because of the piezoelectric effect. The most general so- 
lution of Eq. H3U|) , satisfying the required boundary con- 
ditions, is 

u{z, t) = -b^^ 2 + Re ^ A„ cos(fc„z) e-^"*--"*, (31) 



C33 



where 



fc„ = mr/b. 



(32) 



Here we have assumed that cr is quasi-stationary, so that 
dfa is negligible. The first term in Eq. H31|l describes a 
background strain caused the electric field in the capac- 
itor, present in the classical limit even at zero tempera- 
ture, while the second term describes harmonic fluctua- 
tions about that strain. The n = mode is a center-of- 
mass translation. The n = 1 mode is the fundamental 
thickness-oscillation mode of interest here; it has an an- 
gular frequency given by Eq. (|23|l . 



V 



dz Ez ~ h h^^bU, 

£33 



(33) 



across the resonator and J J, in terms of the charge on 
the electrodes and the spatially averaged strain 



U{t) 



u{b, t) - u{0,t) 



(34) 



in the resonator. Eq. H33|) can then be written in terms 
of the total charge q = (jttR^ on the upper plate as 



q = a,s{V~bh33 U), 



(35) 



where Cros = e337ri?^/6 is the geometric capacitance of 
the resonator. 

The resonator therefore produces a current equal to 



Ircs^aos{V-bh33 U). 



(36) 



The first term in Eq. ((36|l describes a purely capacitive 
effect, which would be present even in the absence of 
the piezoelectric disk between the electrodes. We will 
find that this term simply adds the capacitance of the 
resonator in parallel with the junction capacitance C, 
thereby reducing the junction's charging energy. The sec- 
ond term is a consequence of piezoelectricity, and will be 
shown to have two effects: coupling the JJ to resonator 
phonons and renormalizing Crcs- 

It will be convenient to write Eq. (|31|l as 



m(z, t) 



C33 



z + Su{z, t), 



where 



Su{z, t) = Re \^ An 



cos{mrz/b) e 



(37) 



(38) 



C. Model Hamiltonian 

Next we derive a model Hamiltonian for a single 
current-biased JJ coupled to the dilatational mode of a 
piezoelectric nanomechanical disk resonator. The layout 
is similar to that illustrated in Fig. 0] except that there 
is only one junction, and the gate electrode is not split. 
Extension to multiple junctions and resonators will be 
carried out in Sec. IVII As before, we will assume that 
the junction and resonator states are long lived, and any 
effects of decoherence are neglected. We will proceed by 



is the harmonic fluctuation contribution. After quantiza- 
tion, this latter part of the displacement field will come 
from phonons. The average strain can be similarly ex- 
panded as 



C33 



where 



suit) 



5u{b,t) - 5u{0,t) 



(39) 



(40) 
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Now, the time derivative of the first term in Eq. (|39|l is 
itself proportional to I^^s, so Eq. H3()|l can be equivalcntly 
written as 



where 



Cri 



1 



■7-7 



(41) 



(42) 



is a piezoelectrically enhanced resonator capacitance, and 
7 is the piezoelectric efficiency defined in Eq. H28|l . In 
contrast with that of Eq. H36() , the second term in Eq. H41() 
describes a pure coupling to resonator phonons. 

Returning to the inset of Fig. ^ we replace Ih with 
-^b + -^ros- In our coupled junction-resonator system, /b 
then refers to the bias current coming from the external 
circuitry alone, which may have both DC and RF com- 
ponents (see Sec. IIII A|) . The semiclassical equation of 
motion replacing Eq. ((HJ is now that of a particle with a 
modified mass moving in a potential U + SHd, where 



2e(l — 7 — 7"^) 



(43) 



The classical junction-resonator interaction Hamiltonian 
SHci is evidently linear in the phase difference S. The 
effective mass M of the particle is given by Eq. {Tj) , with 
Ec now reduced to 2e^/(C -I- Cics)- 

Quantization of the S variable proceeds as in Sec. Ill Bl 
The quantization of the resonator dynamics is carried out 
in Appendix^ The resonator Hamiltonian (dropping an 
irrelevant additive constant) is 



(44) 



where and a are bosonic creation and annihilation op- 
erators for dilatational phonons. The junction-resonator 
interaction Hamiltonian is found to be 



SH = —ig{a — a')6 



where 



9 



h^^^ 633 Cr-cs ^/i^ 



(45) 



(46) 



is a real-valued coupling constant with dimensions of en- 
ergy. We note that g depends only on the properties of 
the resonator and is independent of the parameters char- 
acterizing the Josephson junction. The value of g quoted 
in Eq. applies to a fully gated resonator coupled to 
a single J J; for a J J connected to one half of a split- 
gate resonator, such as in Fig. 01 the relevant interaction 
strength is g/2. 

For a fixed disk thickness b, chosen to determine too, 
the interaction strength varies linearly with disk radius 
R. Using the parameters summarized in Table ITTll for a 
15 GHz AIN resonator, we obtain 



g = 2.70 AieV x R[nm], 



(47) 



TABLE IV: Parameters for a single JJ coupled to the res- 
onator of Table IIIII The junction parameters correspond to 
that investigated Ref. 3]. 



critical current Jo 


21 /iA 


Josephson energy Ej 


43.05 meV 


junction capacitance C 


6pF 


charging energy Ec 


53.33 neV 


zero-bias plasma frequency ti;po/27r 


16.4 GHz 


resonant bias current s* 


0.545 


junction-resonator interaction strength g 


0.620 /leV 


resonant vacuum Rabi frequency S7(0) /2n 


8.79 MHz 


resonant Rabi period 27r/Sl(0) 


113.7 ns 



where i?[/im] is the resonator radius in /im. In the sim- 
plest qubit storage simulations carried out below, we 
choose R to be 0.230 ^m, in which case the interaction 
strength is 0.620 fieV. In Table Hvl we summarize this and 
other parameters associated with the most basic coupled 
JJ-resonator system. 

The complete Hamiltonian of the system is 



H = Ho + 5H, with Ho = + H,c 



(48) 



The junction Hamiltonian Hj depends on s, and when s 
is time-dependent, Hj is also time-dependent. We shall 
address this issue below in Sec. IIII Dl Assuming s is 
constant, the stationary states of Hq may be written as 



|mn) = |m),j (g) |n)r 



(49) 



where m — 0,1,2,... labels the junction state and 
n = 0, 1, 2, . . . is the phonon occupation number of the 
resonator. The eigenvalues of Hq are 



hujQ n. 



(50) 



The \mn) and Emn of course depend on s. We will refer 
to the lowest two eigenstates of H^ as the phase qubit, 
and to AE [defined in Eq. H15|l ] as the qubit level spacing, 
even if there are more than two quasibound levels in the 
junction. 

For many applications it is convenient to write the JJ 
Hamiltonian of Eq. I|13() in second-quantized form, as 



Hi 



(51) 



Here cj„ and Cm are creation and annihilation operators 
for the junction states, which can be taken to be either 
fermionic or bosonic because there is only one "particle" 
in the washboard potential. In this same notation, the 
interaction Hamiltonian becomes 



5H 



(52) 



where the Xmm' are dipole-moment matrix elements de- 
fined in Eq. ifT^ll . 
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An important simplification occurs when only the 
qubit states m — 0,1 are included in the JJ. In this case 
the complete Hamiltonian can be written as 



H 




hujoa^a — igifl - 



2^00 a;oi 
,2^10 a;ii, 



I, (53) 



with the JJ operators written as matrices in the {|0), |1)} 
basis. Recall that the diagonal dipole moments Xmm do 
not generally vanish in the current-biased J J, but for s 
not too close to unity are approximately equal to arcsins. 
In the approximation xqq — xh, and dropping an addi- 
tive constant, we can succinctly write (|53|l in terms of 
the Pauli matrices as 

H = —^p-(7z + hjjJao'a — igia — a^)[2:oocro + Xi:,ia,j-], (54) 

where is the identity matrix. Note, however, that 
ccoo 7^ a^ii when s is very close to 1. Finally, when both 
IS.E « fujjQ and g ^ Ai?, the commonly used rotating- 
wave approximation of quantum optics becomes valid. 
Applied to the form (|53|l or (|54() . the Hamiltonian sim- 
plifies to 

i/jc = ^^-Cz + hiuoa'a — igx^i (aa^ — 0^0"+) , (55) 

where <t± = [a^ ± iay)/2. Hjq is the Jaynes-Cummings 
model of quantum optics (written is a basis that is dif- 
ferent from that conventionally used there). 

The Hamiltonian in Eq. H48|) is equivalent to that of 
a few-level atom in an electromagnetic cavity. The JJs 
are analogous to the atoms. The cavity photons here 
are dilatational-mode phonons, which interact electri- 
cally with the junctions via the piezoelectric effect. Cou- 
pling several junctions to a nanomechanical resonator, as 
illustrated in Fig. 13 then makes the system analogous to 
several atoms in an electromagnetic cavity, except that 
here the atomic level spacing and electron-photon inter- 
action strength are all externally controllable. 



D. Quantum dynamics in the instantaneous basis 

As discussed above, the Hamiltonian Hj for the JJ de- 
pends on the dimensionless bias current s, and is there- 
fore usually time-dependent. It will be useful to expand 
the state of the coupled system in a basis of normalized 
instantaneous eigenstates \mn)s of Hq, defined by 

Ho{s)\mTi^^ — Emn{s)\mn'^^, with s = s{t). (56) 

We assume that at time t = to the bias current is constant 
and that the system is prepared in a pure state. For 
t > to we write the wave function, suppressing the time- 
dependence of s{t), as 



^(^)e■ 



Is 



(57) 



The probability amplitudes in the instantaneous interac- 
tion representation satisfy 



{•mn\5H — ihBAm'n' 



Mh)J'dt'lE^„{s)~E^,^,{s)] 



Off-diagonal matrix elements of the quantity 



{mn\-^\m'n') 



(mn\-§^\m'n'% s 



(58) 



(59) 



determine transitions between the instantaneous eigen- 
states caused by nonadiabatic variation of s; the diag- 
onal elements determine the Berry connection of adia- 
batic perturbation theory. In the small s, quadratic- 
potential limit, the low- lying J J eigenstates are well 
approximated by the harmonic oscillator eigenfunctions 
given in Eq. . In this case it can be shown that 



{mn\-^\m'n')_ 



1 ( I m' + l X 

T3if (V^~ 

-l)^nri' 



I 1 rf£W\/(™' + l)(™'+^ S 

- ^^¥^'5™.™'-2)'5™', (60) 



where 



ds 



sL 



4(1 -s2)- 



(61) 



There are no diagonal (Berry connection) terms in this 
limit. The terms in Eq. H60|l proportional to dls/ds re- 
sult from the change of curvature at the minimum of the 
anharmonic potential U{5) with changing s. 



IV. QUBIT STORAGE AND TRANSFER 

We now turn to a discussion of some single-qubit op- 
erations made possible by the nanomechanical resonator. 
In particular, we show that any phase qubit state 



iV-j) =a|0),i+/5|l)j 



with 



a ■ 



1 (62) 



produced in the current-biased JJ can be written to and 
coherently stored in the phonon-number states of the res- 
onator, as 



a|0)re.-)-/3|l) 



(63) 



yielding a quantum memory element. In Eq. (|62() . the 
states |0)j and |l)j are the lowest two junction eigen- 
states shown in Fig. 13 whereas in Eq. (|63|l they denote 
the vacuum and one-phonon states of the resonator's di- 
latational mode. Later, the qubit state can be reversibly 
retrieved or transferred to a second Josephson junction. 

We will examine qubit storage and transfer in two 
stages: First we will develop a simple analytic theory 
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based on the adiabatic approximation combined with 
the rotating-wave approximation (RWA) of quantum 
opticsii^ The adiabatic approximation assumes that the 
bias current s changes slowly on the frequency scale 
AE/h, a requirement that (although not always desir- 
able) can be easily satisfied in practice. The RWA for a 
phase qubit is valid when two conditions are met: 

(i) AE and huiQ are close to each other on the scale 
of the resonator's energy width Hluq/Q. Here Q 
is the resonator's dilatational-mode quality factor. 
Transitions to higher levels with to > 1 are 
far off resonance on this same scale. 

(ii) The interaction strength g is small compared with 
AE (or hujo). 

We will then supplement the analytic theory with numer- 
ical simulations based on the full Hamiltonian of Eq. (|48|l . 
using realistic values of all parameters involved. 



TABLE V: Probability amplitudes Cmn{t) for phase-qubit 
coupled to nanomechanical resonator, at time zero when they 
are brought to perfect resonance, as well as one quarter, one 
half, and three quarters of a vacuum Rabi-oscillation period 
later. 



amplitude 


t = 


t = n/2D. 


t = Tl/fl 


t = 3-K/2D. 


Coo 


Q 


a 


a 


a 


Coi 





(3/V2 


P 


P/V2 


ClO 


13 


P/V2 







Cll 















We have also assumed that all dissipation and deco- 
herence mechanisms are negligible over experimental 
timescales. Furthermore, as discussed in Sec. Ill Bl we 
can take xqi to be positive with no loss of generality. 
Using Eqs. Ht)5|) and (|67|l . we then obtain, by Laplace 
transformation, 



A. RWA analysis 

To understand qubit storage, consider a single junction 
coupled to a nanomechanical resonator as described by 
the Hamiltonian of Eq. (|48|l . and expand the wave func- 
tion for the combined system as in Eq. (|57|l . The proba- 
bility amplitudes Cmn (t) in the instantaneous interaction 
representation then satisfy Eq. (|58() . 

We start at some time Iq < with the JJ prepared in 
the state (I62|l and the resonator in its ground state, 

mo)) = (a|0)j-K/3|l)j)®|0),.es = a|00) + /3|10). (64) 

We assume that the qubit and resonator are detuned and 
that g <C AE. Then the \mn) in (|64(l are close to eigen- 
states, and the Cmn remain approximately constant. The 
qubit level spacing AE is now adiabatically changed to 
the resonant value, reaching hcoo at time t — 0. Then at 
t — we have, approximately. 



Cm„(0) = {a5mQ+ f^S,nl)SnO- 



(65) 



As we shall discuss below, the first nonadiabatic correc- 
tions to Eq. H65|l principally affect the phases of the |00) 
and 1 10) components of the wave function, that is, the 
phases of a and p. The wave function at t = is there- 
fore 



\m) 



^^-(i/K)f,ldtEools(t)] 
P^-{^/ri)f°dtEiols(t)] 



(66) 



Next we invoke the RWA, which allows us to write 
Eq. ^ as 



COn — T 

n 

I 

Cln -- 



^ xoi e*'^'^*ci „_i 



Vn+l xoi e '"■^*co,n-Hi- 



(67) 



coo(i) = a 

co,{t) = /3Msin(f)e--*/2 



cos(f)-)-z^sin(^) 



cioit) = /3 

cii(0 = 0, (68) 
and all Cmn{t) with n > 1 equal to zero. Here 

n{u;d) = ^J[n{0)]^+ujl with r!(0) = ^^, (69) 

is the vacuum Rabi frequency, and Wd = ^^o ~ AE/h is 
the resonator-qubit detuning. fl(0) is the Rabi frequency 
on resonance. Probability amplitudes at selected times 
are summarized in Tabled The wave function at later 
times t > 0, when the system is on resonance, is therefore 

|V(t)) ~ e-'^»[^*l*/''(^ae-(*/'"')'''*>^™[^(*)l|00) 

+ /5sin(§*)e-('/'')'''*>^°^[^(*^l e''^"' |0l) 

+ /3cos(f )e-(*/''^K''*^^°[^(*)le-^"°*|lO)Y (70) 



where s* is the resonant value of the dimensionlcss bias 
current. We emphasize that the result in Eq. (|70|l is only 
approximate. 

After a pulse duration At, the JJ is again detuned from 
the resonator. The final wave functions, in the instanta- 
neous interaction representation, for several important 
choices of At is summarized in Table I VII In this rep- 
resentation the phase factors exp[—{i/h)JdtEm,nis)] ap- 
pearing in Eq. H57|l are suppressed. 

When U At = tt/2, the system is held in resonance 
for one quarter of the vacuum Rabi period, and the final 
state is entangled. In particular, when the qubit begins 
completely in the excited state, a = and /3 = 1, the 
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resulting state is the maximally entangled Bell state 

»J^«|0)„., (71) 

with the resonator returned to the ground state. Simi- 
larly, after three quarters of a Rabi period, the final state 

is 2-5(|01)j- |10)j)(^ |0)ros. 

After half a Rabi period, or At = tt, the phase qubit 
and resonator states are evidently swapped. The cavity- 
QED analog of this operation has been demonstrated ex- 
perimentally in Ref. [t^ . This operation is extremely 
useful in our architecture. In addition to allowing the 
resonator to be used as a quantum memory element, it 
can be used as a quantum bus to transfer a qubit state 
from one JJ to another attached to the same resonator. 

To retrieve a state that has been stored in the res- 
onator, the junction is again tuned to hujQ, except this 
time for one and a half Rabi periods, or H. At — Stt. 
This longer pulse length requirement follows from the 
RWA equations (|67() . An exception is the special case 
where the stored state is a |0) or |1), and one does not 
care about the overall phase of the final result, as in our 
Ref. [sj. We have extensively simulated the use of the 
resonator as a quantum memory element, and its speed 
and fidelity as a function of coupling strength and po- 
sition of the stored state on the Bloch sphere. With 
dimensionless interaction strengths g/hujQ around a few 
percent, the RWA is quite reliable, and the main source of 
error comes from nonadiabatic effects during the ramp- 
ing of s{t), which leads to errors in Eq. I|65(l . These are 
principally phase errors in cqo and cio, which results in a 
significant qubit-state dependence to the memory fidelity, 
with states closer to the equator of the Bloch sphere be- 
ing stored less accurately. We shall return to these issues 
in future work. 

To transfer a qubit state a|0) ji -I- /?|1) ji from junction 
1 to junction 2, the state is stored in the resonator's di- 
latational phonon number states as a|0)ics + /?|l)rcs- Af- 
ter junction 1 is taken out of resonance, the bias on the 
junction 2 is varied to bring it into resonance with the 
resonator for one and a half Rabi periods {fl At = Stt) , 
resulting in the creation of the state a|0)j2-f /3|l)j2 in the 
second junction. (Again, the case where the transferred 
state is a |0) or |1) is exceptional, and a half of a Rabi 
period is sufficient.) The original qubit state is therefore 
transferred from one junction to another. It will be pos- 
sible to verify experimentally that this has occurred by 
reading out the second junction at the end of the transfer 
operation. 

B. Simulating storage and transfer 

The analysis above, which is based on the adiabatic 
and rotating- wave approximations, implies that J J states 
can be stored, transferred, and controUably entangled 
with perfect accuracy, and — with an appropriate choice 
of g — arbitrarily quickly. This is not the case: The actual 



TABLE VI: Approximate final wave functions, in the instan- 
taneous interaction representation, after the phase qubit and 
resonator have been in resonance for a time At. 



nAt 


final state 


operation 


7r/2 


q|00) -F/3(|01) -h|10))/V2 


entangle 


TT 


|0)j(8)(Q|0),es +/3|l)ros) 


swap 


37r/2 


q|00) +/3(|01) - !l0))/\/2 


entangle 



fidelity is determined by the corrections to these approx- 
imations. In this section we shall study the storage and 
transfer fidelities by direct numerical integration of the 
time-dependent Schrodinger equation. 

We begin by simulating the storage of a JJ state in the 
phonon-number states of a resonator. To do this we solve 
the time-dependent Schrodinger equation for the coupled 
junction-resonator system by numerically integrating the 
coupled equations for the case [see Eq. 

a = and (3=1. (72) 

This corresponds to the phase qubit starting in the ex- 
cited eigenstate |l)j. The resonator starts out in its 
ground state |0)res- Our main result, which is shown 
in Fig. will be discussed in detail below. 

To ensure the reliability of the numerical results we 
employed a variety of ODE integrators, including both 
explicit and implicit algorithms, as well as exact diago- 
nalization for cases with constant s. No significant dif- 
ferences were observed. The results presented were ob- 
tained with the 4th-order Runge-Kutta method with a 
time step of 1 fs, which guaranteed that probability was 
conserved for the duration of the calculation to better 
than 99.993%. Josephson junction energy levels and 
dipole-moment matrix elements Xmm' as a function of 
s were calculated using the diagonalization method dis- 
cussed in Sec. IIIBI and found to be extremely close to 
that of a harmonic oscillator in the range of bias currents 
employed here. 

We simulate a large area, current-biased JJ with pa- 
rameters corresponding to that investigated in Ref. [3], 
namely Ej = 43.05 meV and Ec = 53.33 neV. The zero- 
bias plasma frequency u!po/2TT is therefore 16.4 GHz. A 
15 GHz resonator will be in resonance with this junction 
when s = 0.545, comfortably far from the regime near 
s = 1 where bias-current fluctuations are most destruc- 
tive. The nanomechanical resonator we simulate has the 
parameters listed Table UTTl which results in a junction- 
resonator interaction strength g given in Table Hvl The 
resonator thickness d is determined by the desired 15 GHz 
frequency of the thickness-oscillation mode, and the disk 
radius R can be used to vary g without appreciably af- 
fecting that frequency. As we noted in Eq. l|T7|l . g is lin- 
early proportional to R (in the large R/d limit). We have 
used this tunability to ensure that the system is in the 
regime where the RWA analysis of Sec. IIV Al is applica- 
ble. Below we will briefly examine results of simulations 
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TABLE VII: Final state amplitudes Cmn after qubit storage. 
System parameters are the same as in Fig. |S] 



probability amplitude 


Cmn 


Im Cm,i 


j C-mn 1 


Coo 


-0.046 


0.016 


0.002 


coi 


-0.061 


0.992 


0.987 


ClO 


0.049 


-0.007 


0.002 


Cll 


0.045 


-0.030 


0.003 



with larger values of g. There are more than 400 quasi- 
bound states in the junction when s — 0.545. To 
the accuracy of the numerical results reported here, we 
find no sensitivity to the number of JJ states included 
in the calculations as long as at least 4 states are in- 
cluded. The resonator, of course, has an infinite num- 
ber of phonon-number eigenstates |n)rcs, and the results 
shown here have been calculated by including the 4 states 
lowest in energy, as increasing beyond this number led to 
no significant changes. 

We turn now to a discussion of Fig. |S| At time zero 
the current bias is s = 0.40 and the wave function am- 
plitudes are taken to be Cm„(0) = (5mi^nO- The bias is 
held at s = 0.40 for 5 ns. As shown in Fig. |H1 the oc- 
cupation probability of the |10) state remains close to 
unity during this time interval. All other states remain 
essentially unoccupied. After 5 ns the bias current is adi- 
abatically changed to the resonant value of s = 0.545. 
Our simulations show that the success of a qubit storage 
depends somewhat sensitively on the shape of the bias- 
current profile s{t) in the transition region. In particular, 
we find that the time during which s changes from the 
off-resonant value to the resonant one should be at least 
exponentially localized. The result presented in Fig. |^ 
was obtained using a trapezoidal profile with a cross- 
over time of 1 ns, which should be compared with the 
resonator and on-resonance qubit period of 0.1ns. Sim- 
ilar results were obtained using Gaussian profiles. The 
JJ level spacing is tuned to hw^ for half of a Rabi pe- 
riod 7r/f2. During this time interval the junction inter- 
acts strongly with the resonator, and energy is exchanged 
back and forth between the two systems. The JJ is then 
detuned from the resonator. Some of the final probability 
amplitudes are given in Table IVlIl For the small value 
of g used here, chosen so that g/hwo = 0.01, the numeri- 
cal results for the IcmnP are in excellent agreement with 
the RWA. However, the RWA prediction for the phases 
of the Cmn are poor until one goes to even smaller values 
of g. In other words, the RWA is better at predicting the 
moduli of the than their phases. 

It is interesting to examine the extent to which higher 
energy states of the junction and resonator become ex- 
cited during the storage operation. In Fig. Owe plot the 
occupation probabilities of the states |20) and |21), both 
of which involve the higher lying m = 2 junction state. 
Similarly, in Fig. |Slwe plot the occupations of |02) and 



1 12), which involve the n = 2 phonon state. In all cases 
the excitation of higher lying states is negligible. 

A few comments about these results are in order: The 
observed sensitivity to the shape of s{t) can be under- 
stood by recalling that in the absence of any dissipation 
or decoherence, the RWA requires the qubit to be ex- 
actly in resonance with the nanomechanical resonator. 
Therefore it is necessary to bring the two systems into 
resonance as quickly as possible without violating adia- 
baticity. The power-law tails associated with an arctan- 
gent function, for example, lead to considerable devia- 
tions from the desired RWA behavior, as we demonstrate 
in Fig. 1^ We expect this sensitivity to be present in 
real systems as well. We also found that the validity 
of the RWA requires g to be considerably smaller than 
hioo- The ratio g/fuvo for the system simulated in Fig.|2| 
is 1%. When the resonator disk radius R is increased 
to 2.3 /im, g/hujQ is then only 10%, but the RWA al- 
ready fails considerably. This strong-coupling breakdown 
is demonstrated in Fig. ^| The resonant Rabi period in 
this case is 11.4 ns. Of course, the value of At used in 
FigS-l^andlTUIare consequences of the RWA analysis, and 
better fidelity could be obtained by choosing At differ- 
ently. 

Up to this point we have only discussed storage of the 
simple qubit state |1). Storing general qubit states of 
the form a|0) -I- f3\l) follows similarly, although achiev- 
ing high fidelity requires more care. The reason is that 
the ramping of up of s(t) introduces phase errors into 
Eq. the "inital" amplitudes that get swapped. This 
can be circumvented to a considerable extent by choos- 
ing an optimum value of the off-resonant bias current. In 
Fig. ^] we show results of the successful storage of the 
qubit states 2^3 (|0) + |1)) and 2-^(|0) -|-i|l)), which are 
on the equator of the Bloch sphere, using s ~ 0.180 when 
detuned from the resonator. 

Finally, in Fig. 1121 we present results of simulations 
of two junctions coupled to a resonator. The JJs are the 
same as in Fig. El but the resonator in this case has radius 
R = 0.459 /im. Because the upper gate is now split, g = 
0.620 /ieV for each JJ. The instantaneous eigenstates of 
the uncoupled system can be written as \m1m2n), where 
mi and m2 are the eigenstates of the junctions and n is 
the phonon number of the resonator. The phase qubit is 
first stored in the resonator, as described above, and is 
then passed to the second identical junction. The result 
is a transfer of the qubit state |1) from one JJ to another. 
Only half a Rabi period of resonance with the second JJ 
is needed for this transferred state; in general, one and 
a half periods are required. The probability amplitudes 
after the transfer are given in 

Table 



V. TWO-JUNCTION ENTANGLEMENT 

The nanomechanical resonator can also be used to pro- 
duce states where the JJs are entangled, but the res- 
onator remains in its ground state, unentangled with 
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TABLE VIII: Final probability amplitudes Cmj,„2„ after 
transfering qubit state from one junction to another through 
the nanomechanical resonator. Transfer succeeds with a fi- 
dehty squared of better than 97%. 



probability amplitude 


Re C-rni m2 


Im Cm I m2 ^^ 


1 Cm im2'n\ 


cooi 


-0.075 


0.003 


0.006 


coio 


0.591 


0.790 


0.974 


ClOO 


0.023 


0.038 


0.002 



the junctions. We assume that two identical JJs are 
attached to the same split-gate resonator. The instan- 
taneous eigenstates of the uncoupled system are written 
as \m1m2n), where mi and m2 are the eigenstates of the 
junctions and n is the phonon number of the resonator. 

The foundations for this operation have already been 
explained in Sec. IIVI According to Table IVll we can 
prepare an entangled state of two JJs by bringing the 
first junction, previously prepared in the state |l)ji, into 
resonance with the resonator for one quarter of a vac- 
uum Rabi period, or flAt — 7r/2, which produces the 
interaction-representation state 2^2(|001) -I- |100)). The 
first JJ is now maximally enangled with the resonator, 
while the second junction is in the ground state. Af- 
ter bringing the second junction into resonance for half 
of a Rabi period, the state of the resonator and second 
junction are swapped, thereby "passing" the resonator's 
component of the entangled state to the second junction. 
After detuning the second junction, the system is then 
left in the interaction-representation state 



|100)-|010) _ |10)j-|01)j 



V2 



V2 



<»|0)rc 



(73) 



The two Josephson junctions have been prepared in the 
maximally entangled Bell state 2~2(|10)j — |01)j). To 
produce the state 2~2 (|10)j -I- |01)j), the 17 Ai = tt swap 
pulse should be replaced with a 17 At = 37r swap pulse. 

In Fig. El we present the results of a simulation of 
entangled state preparation. The JJs are the same as 
in Fig. IHl and the resonator has radius R — 0.459 /im, 
resulting in an interaction strength g — 0.620 /ieV for 
each JJ. The desired entangled state is prepared with a 
squared fidelity of about 92%. 



couples each adjacent resonator "hub," allowing a com- 
pletely scalable geometry without intrinsic size limits. 

The Hamiltonian for an arbitrary large-scale quan- 
tum information processing circuit consisting of M phase 
qubits and A/" nanomechanical resonators is constructed 
as follows. Let I =1,2, . . . ,Af label the resonators, which 
for simplicity we assume to lie in a two-dimensional plane, 
and let J = 1, 2, . . . , label the junctions. Typically 
there will be many more JJs than resonators. Each junc- 
tion can couple to one or two resonators, subject to the 
constraint that a resonator can support on the order of 10 
junctions, and that, for convenience, bus qubits should 
connect adjacent resonators. The Hamiltonian for such a 
quantum computer, ignoring state preparation, manipu- 
lation, and readout circuitry, as well as all environmental 
coupling, energy relaxation, and decoherence, is 

I Jrn 

- ^ X! X! ~ ^J-mm' ^Im^Jm'- (^4) 

J"-/ mm' 

Here coi is the dilatational mode frequency of resonator /, 
a\ and Oj are dilatational-mode phonon creation and an- 
nihilation operators satisfying [aj,a|,] = Sjji, €jm is the 
spectrum of phase qubit J, and Cj^ and Cjm are creation 
and annihilation operators (either bosonic or fermionic) 
for states m in junction J. The matrix gjj gives the in- 
teraction strength between resonator / and junction J; 
bus junctions have nonzero gu for two values of /, com- 
putational junctions will have only one nonzero element. 
In Eq. (|74|l we have also neglected a small capacitive in- 
teraction between phase qubits connected to the same 
resonator. 

As we will demonstrate in future work, the resonator 
can be used to mediate two-qubit quantum logic between 
phase qubits connected to that resonator. The quantum 
circuit of Fig. 1141 then allows quantum logic to be per- 
formed between any pair of computational qubits Ji and 
J2. This is accomplished by swapping the state stored in 
J2 with a phase qubit J( attached to the same resonator 
as Ji , performing the logical operation on Ji and J[ , and 
then reswapping J( and J2. Any pair of computational 
qubits in Fig. E| can also be controUably entangled. 



VI. LARGE-SCALE QUANTUM CIRCUIT 

A strength of our architecture is scalability: By in- 
troducing additional 6ms junctions coupled to a pair of 
resonators, each resonator with a slightly different di- 
latational mode frequency, the quantum states of the 
resonators can be swapped. This makes it possible to 
construct a large JJ array, with all phase qubits coupled. 
We call this layout a "hub-and-spoke" network, an exam- 
ple of which is shown in Fig. ^1 Each bus qubit "spoke" 



VII. DISCUSSION 

We have introduced a design for a scalable, solid-state 
quantum computing architecture based on the integra- 
tion of nanoclcctromechanical resonators with Josephson 
junction phase qubits. Quantum states prepared in a 
Josephson junction can be passed to the nanomechani- 
cal resonator and stored there, and then can be passed 
back to the original junction or transferred to another 
with high fidelity. The resonator can also be used to 
produce entangled states between a pair of Josephson 
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junctions. Universal two-qubit quantum logic will be ad- 
dressed in future work. The architecture is analogous to 
one or more few-level atoms in an electromagnetic cavity, 
and the junction-resonator complexes can assembled in a 
hub-and-spoke layout, resulting in a large-scale quantum 
circuit. 

The calculations presented here have ignored all effects 
of dissipation and decoherence, with the assumption that 
the associated lifetimes are longer than a few hundred ns. 
This is not unreasonable given the current experimental 
situation. Nor have we attempted to perform the op- 
erations as fast as possible, and we expect there to be 
considerable room for improvement in both speed and 
fidelity. 

Finally, we emphasize that many of our results will ap- 
ply to other resonator- or oscillator-based qubit coupling 

15.17.18.19.20.21.22.23.24.25.26.27.28.29.30.31.32 j^^ 

ticular, the promising design being developed at 
Yaigji5j28j29 ygjjjg charge qubits coupled to superconduct- 
ing transmission line resonators, is very similar to the 
architecture discussed here. 



VIII. ACKNOWLEDGMENTS 

It is a pleasure to thank Steve Lewis, Kelly Patton, 
Emily Pritchett, and Andrew Sornborger for useful dis- 
cussions. MRG was supported by the National Science 
Foundation under CAREER Grant No. DMR-0093217. 
ANC was supported by the DARPA/DMEA Center for 
Nanoscience Innovation for Defence. 



APPENDIX A: QUANTUM MECHANICS OF 
THE PIEZOELECTRIC RESONATOR 

Here we quantize the vibrational dynamics of the piezo- 
electric resonator. In the quantum limit, the first term 
in Eq. (|37|l . which describes the background strain gen- 
erated by the charge a{t), becomes trivially quantized: 
It gets multiplied by the identity operator. 

The quantization of the fluctuation term Su{z, t) pro- 
ceeds similarly to that of ordinary phonons, although we 
have to treat the zero-frequency (n = 0) mode separately. 
First we construct a complete set of orthonormal eigen- 
functions from Eq. 1)38(1 . namely 



/„(z) = -^^-^ cos(n7rz/6), n = 0,l,2,.... (Al) 

These eigenfunctions can be shown to satisfy orthonor- 
mality 

/ dzr:,{z)u^)^5^n (A2) 

Jo 

and completeness 

Y^rjx)Ux')^5(x-x'), (A3) 

n=0 

although in our case the fn{z) are purely real. 



The quantized displacement-fluctuation field is given by 



5u{z) = Mz) zo + J2 y (/»(^) + fni^) 4) , (A4) 

and its associated momentum density 11 = pn^dtu is 

n(z) = fo{z) Po - * E V^^^(/"(^) «" - fni^) 4) , (A5) 

n— 1 ^ ' 

where fc„ is defined in Eq. H32|l. Here zq is the z component of the resonator center-of-mass coordinate operator, 
Po is the z component of the center-of-mass momentum operator, and [^OiPo] = The n = term is excluded in 
the summations of Eqs. (|A4|) and (|A5|) because the corresponding frequency vkn vanishes; its separate inclusion in 
the form given above will enable the use of the completeness relation (|A3() in the analysis below. The and aj^ 
are bosonic phonon annihilation and creation operators satisfying [a„,a|j,] = (S„„'. pun = M-^^s/b is the linear mass 
density of the cylindrical resonator, with M^^s the resonator's mass. Using Eq. (|A3|I it can be shown that 

n(z')] = [<5u(z), n(z')] = ih6{z - z'), (A6) 

as required. 

The final expression for the quantized displacement field is therefore 



u{z, t) = , + ^o(t) + g W— ^ (Uz) a„ e-»'="* + f:{z) al e^^^A , (A7) 

C33 V ^PlinVkn V J 
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where ZQ{t) is in the Heisenberg representation. If the Hamiltonian for the center-of-mass dynamics is Pq/2M^(.s, then 
Zoit) — zq + {po/M^cs)t- Note that the center-of-mass mode does not produce any strain and does not enter into our 
final results. , 



Using Eq. HA7|) leads to 




If we include only the fundamental dilatational (n = 1) 
mode in Eq. (|A8I) . we obtain (suppressing the subscript 
on the dilatational phonon operators) 

SU^-^(a + a^), (A9) 

where i^es = \/fi-/Mrcs is the characteristic size of 
quantum fluctuations in this mode, and where loq is the 
dilatational frequency defined in Eq. H23|) . 

Assuming a harmonic vibrational dynamics for the res- 
onator, and ignoring the center-of-mass motion, the res- 



onator Hamiltonian is 

oo 

H,^s^^hvkn{alan + ^). (AlO) 

n=l 

Keeping only the n = 1 dilatational mode, and dropping 
the additive c- number constant, leads to Eq. (|44|l . Using 
Eq. PTO|) . we then obtain 

SU^iiH^JUi^^Z^l^i'^-'iy (All) 

n odd 

The n — 1 term on the right-hand-side of Eq. ljAll|l . 
when inserted into Eq. (|43|l , yields the interaction Hamil- 
tonian of Eq. (|45|l with the coupling constant given in 
Eq. 
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FIG. 1; (main panel) Effective potential U{S) for dimensionless bias current s = h/Io equal to 0.1, plotted in units of Ej. 
(inset) Equivalent-circuit model for a current-biased Josephson junction. A capacitance C and resistance R are in parallel with 
an "ideal" Josephson element, represented by a cross and having critical current Jq. A bias current lb is driven through the 
circuit . 

FIG. 2: Barrier height and plasma frequency as a function of the dimensionless bias current s. Here AUq = 2E] is the barrier 
height at zero bias, and tJpo is the zero-bias plasma frequency defined in Eq. 1121 . 



FIG. 3: Metastable potential well in the cubic limit, showing the barrier of height AU that separates the metastable states 
|0), |1), and 1 2), from the continuum. This figure applies to the case of bias currents s just below 1. The lowest two states are 
separated in energy by AE. 
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FIG. 4: Two-qubit circuit diagram. The computational qubits are the two JJs in the center, shown as crossed boxes, each 
coupled to one side of the piezoelectric disk resonator. Each crossed box represents a real JJ, modeled by an ideal Josephson 
element in parallel with a resistor and capacitor. The current bias and readout circuits for each qubit circuit are shown on the 
left and right sides of the figure. Note that there is no direct electrical connection between the two qubits. 



FIG. 5: Four current-biased JJs coupled to a nanoelectromechanical resonator. Each junction is connected to a metallic 
plate on the surface of the resonator that covers about one quarter of the surface. Because we make use of the fundamental 
dilatational mode, which is spatially uniform in the plane of the resonator, the qubits are all equally well coupled to that mode. 



FIG. 6: Phase qubit storage. The solid descending curve is |cioft)P, the interaction-representation occupation probability of 
the 1 10) state, calculated numerically for the junction of Ref. |3| coupled to the 15 GHz piezoelectric resonator described in 
Table UTTI The dashed curve is the same quantity calculated from the analytic RWA results of Sec. IIV K\ The solid ascending 
curve is |coi(i)p. The dotted curve shows the time dependence of the dimensionless bias current s{t), which is varied to bring 
the phase qubit in resonance with the resonator after 5 ns. The Rabi period on resonance, when s — 0.546, is 113.69 ns. After 
the storage operation, |ciop = 0.002 and jcoi]^ — 0.987. The inset shows an enlarged view of |cio(t)P during the ramping up 
of s(t). 



FIG. 7: Occupation of higher lying m — 2 junction state during qubit storage. The upper plot is jc2oP, and the lower plot is 
jc2ij^. Both quantities would vanish in the RWA. All junction and resonator parameters are the same as in Fig. |S| 



FIG. 8: Occupation of higher lying n = 2 resonator state during qubit storage. The upper plot is jco2p, and the lower plot is 
|ci2p; both vanish in the RWA. Parameters are the same as in Fig.|S| 



FIG. 9; Qubit storage with arctangent bias-current profile. All system parameters are the same as in Fig. |S| The numerical 
result for |cioP, shown as a solid descending curve, is entirely different than that predicted by the RWA (dashed curve), even 
though the difference between the s{t) profiles used here and in Fig. |S|is small. The qubit state is not correctly stored in the 
resonator. 



FIG. 10: Qubit storage in larger resonator. Here we simulate qubit storage in a 15 GHz resonator with 7? = 2.3 /xm, so that 
g/fkoo = 0.10. All other resonator and junction parameters are the same as in Fig. |n] The solid descending curve is |ciop, and 
the dashed curve shows the desired RWA behavior. The solid ascending curve is |coip. The RWA breaks down here because of 
the stronger interaction strength. The dotted curve is s{t). Qubit storage fails again. 

FIG. 11: Storage of qubit states on the equator of the Bloch sphere, (a) Here the initial state is 2~ 2 (|0)j + 1 l)j) (g) jO)rcs. 
The solid descending curve is the squared overlap with the interaction-representation state 2^2 (|00) -|- \W)), and the ascending 
curve is the occupation of 2^2 (|00) -I- |01)). The dotted curve is s(t). (b) The initial state is 2~ 2 (|0)j +z|l)j) ® |0)rcs. The 
descending and ascending curves are the occupations of 2~2 (|00) -|- i|10)) and 2~2 (|00) + j|01)), respectively. 

FIG. 12: Qubit transfer between two identical Josephson junctions. The descending solid curve is |cioo(t)P, the probability 
for the first junction to be in the m = 1 excited state, and the rest of the system to be in the ground state. The state of the 
first junction is stored in the resonator as in Fig. |S1 the peaked curve giving |cooi(i)p. The ascending curve is |coio(t)|^, the 
probability for the second J J to be in the excited state. The solid and dotted trapezoidal curves show the bias currents Si(t) 
and S2 (t) on the two junctions. 

FIG. 13: Preparation of entangled Josephson junctions. The thick solid curve is the probability for the system to be found in 
the interaction-representation state 2~2 (|100) — |010)). The thin solid and dashed lines are Si{t) and S2{t), respectively. 

FIG. 14: Architecture for a large-scale JJ quantum computer. In addition to the junctions coupled to a single resonator, as in 
Fig. |S] here there are additional bus junctions for transferring states between different resonators. 



This figure "figl.jpg" is available in "jpg" format from: 



http ://arXiv . org/ps/quant-ph/0409 1 7 9 v 1 



This figure "fig2.jpg" is available in "jpg" format from: 



http ://arXiv . org/ps/quant-ph/0409 1 7 9 v 1 



This figure "fig3.jpg" is available in "jpg" format from: 



http ://arXiv . org/ps/quant-ph/0409 1 7 9 v 1 



This figure "fig4.jpg" is available in "jpg" format from: 



http ://arXiv . org/ps/quant-ph/0409 1 7 9 v 1 



This figure "fig5.jpg" is available in "jpg" format from: 



http ://arXiv . org/ps/quant-ph/0409 1 7 9 v 1 



This figure "fig6.jpg" is available in "jpg" format from: 



http ://arXiv . org/ps/quant-ph/0409 1 7 9 v 1 



This figure "fig7.jpg" is available in "jpg" format from: 



http ://arXiv . org/ps/quant-ph/0409 1 7 9 v 1 



This figure "fig8.jpg" is available in "jpg" format from: 



http ://arXiv . org/ps/quant-ph/0409 1 7 9 v 1 



This figure "fig9.jpg" is available in "jpg" format from: 



http ://arXiv . org/ps/quant-ph/0409 1 7 9 v 1 



This figure "figl0.jpg" is available in "jpg" format from: 



http ://arXiv . org/ps/quant-ph/0409 1 7 9 v 1 



This figure "figll.jpg" is available in "jpg" format from: 



http ://arXiv . org/ps/quant-ph/0409 1 7 9 v 1 



This figure "figl2.jpg" is available in "jpg" format from: 



http ://arXiv . org/ps/quant-ph/0409 1 7 9 v 1 



This figure "figl3.jpg" is available in "jpg" format from: 



http ://arXiv . org/ps/quant-ph/0409 1 7 9 v 1 



This figure "figl4.jpg" is available in "jpg" format from: 



http ://arXiv . org/ps/quant-ph/0409 1 7 9 v 1 



